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ABSTRACT 


On  the  basis  of  assumed  additivity  of  risk  functions,  the 
statistical  distribution  cf  the  time  to  first  failure  of  a  group 
of  nominally  identical  structures,  each  with  two  critical  loca¬ 
tions  respectively  subjected  to  chance  failure  and  fatigue  fail¬ 
ure,  is  derived  together  with  analytical  expressions  of  other 
statistical  quantities  pertinent  to  the  problem  of  the  time  to 
first  failure.  The  concept  of  fatigue  sensitivity  is  examined. 

An  example  involving  aircraft  wings  designed  according  to  che 
current  method  of  design  is  worked  out  to  demonstrate  how  the 
reliability  of  such  wings  can  be  estimated.  The  example  shows 
a  dominance  of  chance  failure  if  the  group  size  increases  and/cr 
a  higher  reliability  is  demanded.  It  is  pointed  out  that  the  pre¬ 
sent  method  of  reliability  estimation  can  easily  be  extended  to 
cover  mere  realistic  structuras. 


This  abstract  is  subject  to  special  export  controls  and  each 
transmittal  to  foreign  governments  or  foreign  nationals  may 
be  made  only  with  prior  approval  of  the  Metals  and  Ceramics 
Division  (MAM),  Air  Force  Materials  Laboratory,  Wright- 
Patterson  AFB,  Ohio. 
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SECTION  I 


INTRODUCTION 

Significance  of  the  time  to  lirst  failure  (minimum  life  among 
a  group  of  nominally  identical  structures  each  subjected  to  a  sam¬ 
ple  sequence  of  the  random  load)  as  one  of  most  important  design 
parameters  from  the  view  point  of  structural  reliability  is  empha¬ 
sized  by  A.  M.  Freudenthal. ^  The  emphasis  is  essentially  based  on 
the  following  grounds:  (1)  The  implication  of  a  high  percentage 
of  number  of  structural  failures  before  the  mean  time  to  failure 
precludes  the  use  of  the  mean  time  to  failure  as  a  design  parameter, 

(2)  the  time  to  first  failure  is  of  major  concern  in  the  engineer¬ 
ing  as  well  as  economical  and  other  decision  making.,  (3)  mathemati¬ 
cally  it  has  a  statistical  distribution  of  extremes  and  hence  the 
use  of  its  asymptotic  form  is  justified  if  a  large  number  of  struc¬ 
tures  are  involved,  and  (4)  its  reliability  demonstration  is  not  com¬ 
pletely  out  of  question  even  with  maj^r  structures  such  as  an  aircraft. 

If  the  structure  to  be  considered  has  only  one  critical  location 
with  a  mode  of  failure  identified  by  a  power  law  of  failure  rate 
(chance  failure  is  a  particular  case  of  this) ,  or  with  a  mode  of  fail¬ 
ure  associated  with  the  log-normal  or  the  gamma  life  distribution, then 

the  statistical  distribution  of  the  time  to  first  failure  among  a  group 

2 

of  these  structures  can  be  computed  without  much  difficulty. 

It  is,  however,  evident  that  such  a  simple  model  is  not  valid  for 
a  structure  such  as  an  aircraft  wing  where  a  large  number  of  locations 
structurally  sensitive  to  fatigue  damage  as  well  as  chance 
failure  have  to  be  considered.  Furthermore,  the  effect  of  interaction 
between  such  chance  and  fatigue  failure  on  the  time  to  first  failure  is 
not  clearly  known  at  the  present  time. 

In  the  following  section  (Section  II) ,  therefore,  a  clear  analyti¬ 
cal  formulation  of  the  problem  is  given  and  the  effect  of  interaction 
between  these  two  modes  of  failure  on  the  time  to  first  failure  is  iden- 
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tified.  In  doing  so,  a  structural  model  is  considered  that  con¬ 
sists  of  two  critical  failure  components,  one  in  chance  failure  and 
one  in  fatigue.  This  is  an  idealization  of  the  situation  in  which 
the  root  of  an  aircraft  wing  has  a  possibility  of  failure  in  insta¬ 
bility  on  the  upper  surface  whereas  it  is  sensitive  to  fatigue  fail¬ 
ure  on  the  lower  surface. 

In  Section  III,  the  definition  of  fatigua  sensitivity,  a  notion 
originated  by  A.  M.  Freudenthal^  is  reexamined  and  the  chance-fail¬ 
ure  sensitivity  is  also  introduced. 

Taking  the  aircraft  wing  mentioned  above  as  an  example,  it  is 
demonstrated  in  Section  IV,  (1)  how  the  distribution  of  the  time 
to  first  failure  among  a  group  of  n  such  wings  designed  according 
to  the  current  method  of  design,  can  be  estimated,  (2)  what  effect 
on  the  reliability  is  to  be  expected  if  the  group  size  is  increased 
and  (3)  what  reliability  statements  can  be  made  in  general. 

The  problem  of  how  to  deal  with  more  realistic  cases  involving 
a  number  of  critical  locations  is  discussed,  among  other  items,  in 
Section  V,  suggesting  that  the  method  developed  here  can  be  applied 
to  the  reliability  estimate  of  real  wing  structures  if  these  criti¬ 
cal  locations  and  their  modes  of  failure  as  well  as  the  failure  rates 
are  identified  and  assessed  by  expert  engineers  familiar  with  struc¬ 
tural  details  on  the  basis  of  available  analytical  and  experimental 
knowledge  combined  with  their  experience. 

It  is  emphasized  that  the  reliability  estimate  of  a  structure 
cannot  be  better  than  the  present  state  of  art  of  structural  analysis, 
prediction  of  its  operating  conditions  and  estimation  of  mechanical 
properties  of  the  material  used.  The  fact  that  a  number  of  explicit 
and  quantitative  assumptions  have  to  be  made  in  order  to  arrive  at  a 
final  reliability  figure,  is  not  to  be  considered  as  a  weakness  in¬ 
herent  to  the  reliability  analysis.  Rather,  it  is  to  be  considered 
as  a  clear  and  conscientious  recognition  of  the  problem  areas  in  which 
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further  study,  analytical  and  experimental,  and  statistical  data 
collection  are  needed.  In  fact,  it  is  only  through  the  relia¬ 
bility  analysis  that  the  importance  of  these  assumptions  or  the 
significance  of  additional  information  can  be  assessed  in  direct 
relation  to  the  structural  safety. 
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SECTION  II 


SASIC  PROBABILISTIC  ANALYSIS 


Let  C  |F)  denote  a  chance  failure  (a  fatigue  failure)  at 

location  0  (0  )  within  the  structure.  Furthermore,  let  h  (t) 

C  F  C 

(^h  (t)y  denote  tie  risk  function  or  the  rate  of  failure  for  the 

failure  mode  C  (the  failure  mode  F) .  By  definition, h  (t) h  (t)  '  is 

c  f  y 

such  that  h  (t)dt(h  (t)dt)  is  the  conditional  probability  of  fail- 
C  F 

ure  under  C (F)  in  the  time  interval  [t,t+dt)  given  that  0  (0  )  has 

C  F 

survived  the  interval  [0,t). 


The  chance  failure  implies  a  constant  risk  and  hence 


hc(t) 


in 


whereas  it  is  assumed  as  usual  on  the  basis  of  available  data  that 
the  fatigue  failure  is  a  process  with  the  following  risk  function 
with  a  =  2  ~  5. 


hF(t) 


_a 

v 


L 

vv  J 

F  F 


a-1 


(2) 


where  V  and  V  are  the  characteristic  lives  of  0  and  0  respec- 
C  F  C  F 

tively  and  the  scale  parameter  a  is  an  inverse  measure  of  dis¬ 
persion  of  the  fatigue  life  distribution. 


It  follows  from  the  definition  of  the  risk  function  that 


d  Xn  L  (t) 

d  L  (t) 

hF(t>  ■  ‘  at  * 

_  _ E _ 

dt 

■/LF(t) 

(3) 

d  Xn  L  (t) 

d  L  (t) 

3* 

n 

ri¬ 

ll 

1 

a 

rt 

r 

ii 

C 

dt 

/Lc(t) 

(4) 

„t 

t  “ 

L  (t)  =  exp [-  j  h  (u)du] 

=  exp [~ 

(-)  1 

(5) 

"o 

F 

=  probability  that 

0F  will 

survive  the  interval 

o 

rt 

or,  more  rigorously,  probability  that  a  chance-fail¬ 
ure-free  reference  structure  will  survive  [0,t). 
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(6) 


L  (t)  =  exp  [-  h  (u)  du.]  =  exp  [-(—)] 

C  -'o  C  VC 

=  probability  that  will  survive  the  interval  [0,t) 

(survival  of  a  fatigue-failure-free  reference  structure) 

where  the  initial  condition  L  (0)  =  L  (0)  =  1  is  assumed. 

F  C 

Since  the  survival  of  the  structure  requires  the  survivals 

of  both  0  and  0  , 

C  F 


L (t)  =  Lp(t)  .Lc(t)  =  exp[-  j  (hF(t)+hc(t)  }dt] 


(?) 


=  exp[-  l^~)  -  (“-)  ] 
C  F 


=  probability  of  survival  of  the  structure 
(as  a  whole) in  (0,t) 


From  Eq.  (7),  it  follows  that  the  failure  rate  h(t)  of  the 
structure  is 

h(t)  =  hc(t)  +  hp(t)  =  —  +  f-  (f-)  (S) 

C  F  F 

This  indicates  that  the  probability  of  failure  of  the  struc¬ 
ture  in  the  time  interval  [t,t+dt)  (due  to  the  failure  mode  C  and/ 
or  the  mode  F)  is  the  sum  of  the  probabilities  of  failure  at  O^,  and 
0^  in  the  same  time  interval.  Therefore,  it  implies  that  the  fail¬ 
ures  at  0  and  0  are  mutually  exclusive.  It  is  pointed  out  that 
C  F 

such  mutual  exclusiveness  follows  from  the  independence  assumption 

of  the  failure  processes  at  0  and  0  tacitly  employed  in  deriving 

C  F 

Eq.  (7). 

The  assumption  of  independence  and  hence  the  assumption  of 

the  mutual  exclusiveness  is  not  entirely  correct  since  0  and  0  are 

C  F 

the  parts  of  the  same  structure  and  therefore  are  subjected  to  the 
same  sequence  of  load.  Nevertheless,  Eq.  (8)  is  assumed  in  the  pres¬ 
ent  study  because  (1)  it  is  an  upper  bound  of  the  failure  rate  of  the 
structure  (this  can  be  shown  by  a  direct  use  of  the  addition  rule  of 
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probability)  ani  (2)  it  simplifies  the  analysis  significantly.  The 
following  ana  /sis  provides  a  conservative  reliability  estimate 
because  such  m  upper  bound  is  used  for  the  failure  rate  of  the 
structure . 

The  mo/cality  distributions  f  (t)  and  f  (t)  under  failure 

F  C 

mode  F  and  C  are  respectively  defined  such  that 

a-1  a 

fp(t)dt  =  L(t)h  (t)dt  =  f-  (~ )  exp  [-  (~)  -  (“)  )  (9) 

F  F  C  F 

=  probability  that  the  structure  will  fail  at  0^ 
in  (t,t+dt) 

a 

fclt..  it  =  L(t)h  (t)dt  =  —  exp[-(^-)-  (— )  ]  (10) 

C  C  F 

=  probability  that  the  structure  will  fail  at  C>c 
in  [t,t+dt) 

whereas  tie  mortality  distribution  f(t)  of  the  structure  regardless 
of  the  n.  -de  of  failure  is 

f(t,  =  f  (t)  +  f  (t)  =  L (t) h (t) 

F  C 

=  tr  +f-  <^r,a  1]  e*p>-<r>-  (11) 

C  F  F  C  F 


Define 


P„  =  !  f  (t)dt 
F  v  F 


=  probability  that  the  structure  will  fail  at  O 


pc  ■  J  £c(t)dt 


=  probability  that  the  structure  will  fail  at  0 


Then 


CD 

« 

,  P  +  r  =  f(t)dt  and  therefore 
F  C  J 


P  t  P  =1 
F  C 


This  indicates  the  condition  that  the  failure,  whether  at  0^. 

or  at  0  will  eventually  occur. 

F 


Let  T  denote  a  random  variable  representing  the  life  of  the 
structure.  Then,  the  (conditional)  expectation  of  T  of  the  struc¬ 
ture  which  will  fail  due  to  fatigue  is 

oo 

V  =  j  t  fp(t)dt  /  PF  (15) 

o 

Similarly,  the  (conditional)  expectation  of  T  of  the  structure 
which  will  fail  under  mode  C  is 

CG 

E  T  =  f  t  f  (t)dt  /  P  (16) 

C  J  C  C 

o 

The  expected  value  of  the  life  of  the  structure  is 

CO 

ET  =  f  t  f ( t) dt  =  P  .E  T  +  P  .E  T  (17) 

J  FFCC 

o 

By  performing  integration  by  parts  in  Eq.  (17)  for  the  particular 
form  of  f(t)  given  in  Eq.  (11),  one  can  shew  that 

ET  =  V  ?  (18) 

C  C 

Note  that  V is  the  expected  life  of  a  (hypothetical)  fatigue-free 

"reference"  structure  and  therefore  P  indicates  the  reduction  in 

C 

the  expected  life  (in  fraction)  by  adding  a  location  (or  a  struc¬ 

tural  component)  sensitive  to  fatigue  failure  to  the  reference  struc¬ 
ture.  Probabilistically,  this  is  a  "series  combination"  of  components 
since  the  structure  is  assumed  to  fail  if  one  of  these  components  fails. 

To  introduce  the  concept  of  the  time  to  first  failure,  consider  a 
group  of  n  structures,  each  subjected  to  an  independent  sequence  of 
load,  and  define 

Lp(n)(t)  =l'LF(t)}n  (19) 

=  probability  that  n  chance-failure- 
free  reference  structures  will  all 
survive  the  interval  [0,t)  or  equi¬ 
valently  probability  -that  the  first 
failure  of  these  n  reference  struc¬ 
tures  will  not  occur  in  [0,t) 
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(20) 


Lc(n)(t)  =  {Lc(t)]n 

=  probability  that  n  fatigue-failure-free  structures 
will  all  survive  the  interval  [Ojt)  or  equivalently 
probability  that  the  first  failure  of  these  n  refer¬ 
ence  structures  will  not  occur  in  [0,t) 

L(n)(t)  =  Lp(n)  (t)Lc(n)  (t)  =  [LF(t)Lc(t)]n  (21) 

=  probability  that  the  first  failure ,  whether  under 
mode  F  or  C  will  not  occur  in  [0,t) 


It  follows  from  Eqs.  (19-21)  that  the  risk  functions  h  '  (t) 


and  h  '  (t)  for  the  first  failure  at  0  and  0  respectively  and 

/j.  \  C  f  c 

h'  (t)  for  the  first  failure  regardless  of  the  failure  mode,  are 


given  as  follows. 


(t)  =  nh  (t)  -  - 
F 


(t)  =  nhc (t)  =  - 


d  in  Lp(n)  (t) 


a-1 

na  .  t  . 

V  V 
F  F 


d  in 


T  (n) 


h^  (t)  =  nh  (t)  =  n[h  (t)  +  h  (t)  ] 

r  C 

(n)  (24) 

d  in  Lv  '  (t)  t a-1 

_ _ n_  na  .  t 

dt  ”  V  +  V  lvJ 

C  F  F 

in  which  h  ^(t)dt(h  ^(t)dt)  is  the  conditional  probability  that 
F  C 

the  first  failure  will  occur  in  [t,t+dt)  under  mode  F(0)  given  that 
none  of  n  chance-failure-free  (fatigue-failure-free)  reference  struc¬ 
tures  has  failed  in  the  interval  [0,t)  whereas  h^(t)dt  is  the  con¬ 
ditional  probability  that  the  first  failure  will  occur  in  the  same 
interval  [t,t+dt)  whether  under  mode  F  or  C  given  that  none  of  the  n 
structures  has  failed  in  the  interval  [0,t). 

The  mortality  distributions  f  ^  (t)  and  f  ^  (t)  under  mode  F 

F  C 

and  C  are  defined  such  that 
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(25) 


f  (n)(t:)dt  =  L(n)  (t)h  (n)  (t)dt  =  n{L(t))nhIlt)dt 


-  probability  that  the  first  failure  will  occur  at 

0  in  [t,t+dt) 

F 


fc  (n)(t)ht 


L(n)  (t)hc(n)  (t)dt  =  n{L(t)  }nhc(t)dt 


=  probability  that  the  first  failure  will  occur  at 
0^  in  [tjt+dt) 

whereas  the  mortality  distribution  f ^  (t)  of  the  first  failure 
regardless  of  the  mode  is 

f(n)(t)dt  =[fF(n)  (t)  +  fc(n)(t))dt=  n[L(t))nh(t)dt  ( 2 ' 


Then, 


=  probability  that  the  first  failure  will  occur 
either  at  0^  or  at  0^  in  [t,t+dt) 

p°  (n)  ,  na  p"  t”'1  ~  n  [ ‘'T1  +  (^>  ! 

-  j  fF  <t)dt  =  —  J  (— )  e  c  dt  (28) 

o  F  O  F 


=  probability  that  the  first  failure  will  occur 
under  mode  F 

t  t  Q 

P*  („)  „  r*  -"((•»)  +  (—: l  1 

-  J  *c,n,U)dt  -  J-J  «  *  F  dt 

o  CO 

=  probability  that  the  first  failure  will  occur 
under  mode  C 


As  before. 


(n  n 

p  V  #  p  V  / 

F  C 


=  J  (fF(n)  (t)  +  fc(n)  (t)]dt  =  1 


indicating  that  the  first,  failure  will  occur  eventually. 


Let  Tv  denote  a  random  variable  representing  the  time  to 
first  failure  among  the  group  of  n  structures.  Then,  the  (condi¬ 
tional)  expectation  of  of  the  group  in  which  the  first 
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failure  occurred  under  mode  F  is 


(n)  p  (n)  ,  (n) 

E  T  =  I  t  f  Udt/P  U 

F  J  F  F 

o 


Similarly,  the  (conditional)  expectation  of  T  of  the 

group  in  which  the  first  failure  occurred  under  mode  C  is 


E  T^n‘  =  f  t  f  ^dt  /  P 
C  w  C  C 


The  expected  value  of  the  time  to  first  failure  is 


=  J  t  f(n)(t)dt  =  PF(n.>EFT,n)+  Pc('.')E/,T(n) 


For  the  particular  forms  of  risk  functions  assumed  in  Eqs.  (1) 
and  (2) ,  it  can  be  shown  that 

ET<n)  „  !c  p  (n)  ,  ET  (»)  (n) 

n  C  C  C 

where  the  well-known  fact  that  V^/n  is  the  expected  time  to  first 

failure  in  a  group  of  n  hypothetical  fatigue-free  "reference"  struc- 

(n) 

ture,is  used.  Hence,  Pc  indicates  the  reduction  in  terms  of  frac¬ 
tion  in  the  expected  time  to  first  failure  among  the  group  by  adding 
a  location (or  a  structural  component)  sensitive  to  fatigue  failure  to 
each  of  these  reference  structures. 

It  is  also  pointed  out  that  by  setting  n  -  1,  the  expressions 
given  in  Eq.  (19)-(34)  produce  the  obvious  identities  L>F(t)  =  L  ^  (t)  , 
Lc(t)=  Lc(1)  (t),  hp(t)=  hF<1)(t),  hc(t)=  hc^)  (t),  etc. 

In  the  following  analysis,  the  probability  P  ^  plays  an  impor¬ 
tant  role  and  therefore  is  evaluated  here. 

Setting  f  =  t/V  ,  y  =  V  /V  and  X  =  \n^  in  Eq.  (29),  the 

C  F  C 

probability  can  be  written  as 


=  xj  r‘«  di  =  J  e-1’-<T'/x)  dl 


For  a  =  2,  the  integrals  in  Eq.  (35)  can  be  carried  out  and 
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p  =vtt  \e  [1  -  $  v~  )]  (3o) 

C 

where  $ (x)  is  the  standardized  normal  distribution  function. 

For  a  -  2,  using  Eq.  (36)  and  for  a  =  3,  4  and  5,  evaluating 
the  integral  numerically  on  an  IBM  7094,  the  value  of  ^  is  com 
puted  as  a  function  of  X  and  plotted  in  Fig.  1.  Also,  the  values 
of  P ^  for  various  values  of  n,  a  and  y  are  computed  and  listed 
in  Table  I. 

For  small  values  of  X,  P^n^  can  be  approximated  by  taking 
first  a  few  term  of  the  following  asymptotic  expansion  (take  the 
second  member  of  Eq.  (35),  expand  e  5  and  integrate  term  by  term) 


p  (n)=  --  y  (-D 

C  a  L 


k+1  X  r  .k 

(k-i)  !  *  V 


(1/X  -  ») 


with  T(- )  being  the  gamma  function,  whereas  for  large  values  of  X, 
it  has  the  following  asymptotic  expansion  (take  the  third  member 
of  Eq.  (35),  expand  and  integrate  term  by  term), 

co  k+1 

pc(n)  =  1  -  y  -  r  U+te)  (*  -  »)  (38) 

k=l  k:  X 


Similar  asymptotic  expansions  can  be  obtained  for  E  T  and 
FT1"1. 

ef  ; 

(n )  » 

p  (n)  ■  -c-  ■  7  r  =-  y  -  r(— )  (i/x  -  »)  (39] 

C  (n)  a  L  (k-1) :  'a 

Cj  a  _  ,  - 

C  k=l 


/  ixk+1 

t  _  v  _ 

L  ,  ka 

k=i  k-  x 


f  { 2  +  ka )  ( X  -  00 ) 


E  T(n) 

(n)  V__ 


r-i)^1  xk 

(k-i) : 


f(l  +  k/a)  (1/X  - 
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It  can  be  shown  that  the  compatible  use  (Eqs.  (37) ,  (39)  and 

(41)  together  or  Eqs.  (38),  (40)  and  (42)  together)  of  these  asymp¬ 

totic  expressions  satisfy  the  identity  given  in  Ea.  (33). 


In  Appendix,  the  expression  for  the  coefficient  of  variation 
is  derived  and  computed. 


of 
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SECTION  III 


FATIGUE  SENSITIVITY 


The  concept  of  the  fatigue  sensitive  structure  was  first  in¬ 
troduced  by  Freudenthal^  who  indicated  a  number  of  possible  ways 
of  defining  the  fatigue  sensitivity  factor.  For  example,  take  a 
structure  which  is  subjected  to  two  modes  of  failure  as  discussed 
in  the  preceding  section.  Consider  the  expected  life  ET^  of  this 
structure  disregarding  the  possible  failure  under  the  failure  mode 
F  (hypothetical  fatigue-free  reference  structure)  ; 


ET  =  V  (43) 

C  C 

Similarly,  consider  the  expected  life  ET  of  the  same  structure 

F 

disregarding  the  possible  failure  under  the  mode  C  (hypothetical 
structure  free  of  chance  failure) ; 


ET  =  V  HI  +  1/a) 
r  r 


(44) 


Then,  it  is  possible  to  define  ET^/ET  as  a  measure  of  fatigue 
t  C  F 


sensitivity. 


However,  a  different  definition  of  the  fatigue  sensitivity  can 
be  introduced  through  the  following  ratio  k^. 

kC  =  <£V  ET)/ETC=  1  “  ET/ET„  (0  s  k  <;  1)  (45) 

where  ET  is  the  expected  life  of  the  structure  as  defined  in  Eq.  (17) . 

The  significance  of  this  ratio  is  rather  evident;  shorter  the 
expected  life  of  the  actual  structure  subjected  tu  both  fatigue  and 
chance  failure  compared  with  that  of  the  reference  (fatigue-free) 
structure,  the  larger  k  .  Two  limiting  cases  are  of  interest;  if 
ET  =  0,  then  k^=  1  (completely  sensitive)  and  if  ET  =  ET^,,  then 
k^=  0  (completely  insensitive) . 

Using  Eqs.  (18),  (43)  and  (45), 

kc  =  1  -  pc  (46) 
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which  defines 


Similarly,  one  can  consider  the  ratio  k 


(n) 


the  fatigue  sensitivity  of  a  group  of  n  structure; 

k  (ET  (n)-  ET(n))/ET  (n)  (47) 

v>  V  L 

where  ET^^n^  is,  as  defined  previously,  the  expected  value  of  the 
time  to  first  failure  within  the  group  of  n  reference  (fatigue-free) 
structures . 

Using  Eqs.  (34)  and  (47)  .  artc>  the  fact  that  ET  v  /n, 

kc(n)  =  1  -Pc(n)  ('8) 

Here  again,  an  identity  k^=  k^  ^  is  obtained  (since  P^=  P^1^) 
as  it  should  when  one  set  n  =  1  in  Eq.  (47) . 


The  effect  of  the  group  size  n  on  the  fatigue  sensitivity  k 


(n) 


is  shown  in  Fig.  2  for  y  =  1.0,  0.1  and  0.01  which  indicates  that 
with  all  other  parameters  being  identical,  the  sensitivity  decreases 
as  any  one  of  n,  y  and  a  increases.  However,  when  y  «  1  (as  in  the 
later  example  in  Section  IV),  it  seems  more  convenient  to  consider 
"chance-failure"  sensitivity  k  ^  rather  than  " fatigue" sensitivity 
k^1  since,  with  V  «  1,  the  life  of  the  structure  is  essentially 
controlled  by  fatigue  failure  process.  For  this  purpose,  define 

(49) 


k_^  =  (ET  -  ET^)/ET  ^ 

F  F  F 


where  ET  ^  is  the  expected  time  to  first  failure  among  a  group  of 
F 

n  hypothetical  "chance-failure-free"  reference  structures,  and  it 
can  be  shown'*’  to  be 

(n)  -1/a  _  -1/a 


ET. 


=  n 


ETf=  n 


V  r(l  +  l/a) 
F 


Hence 


k  (n)  =  1  -  PCW/(\T(1  +  1/a)) 


(50) 


(51) 


Equation  (49)  indicates  that  k  ^ET  ^  is  the  reduction  in 

F  F 

the  expected  time  to  first  failure  among  a  group  of  n  chance-fail- 
ure-frea  reference  structures  due  to  an  addition  (in  series)  of  a 
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t 


structural  component  sensitive  t  chance  failure  to  each  of  these 

n  structures.  Therefore  the  closer  k  ^  to  unity*  the  larger  the 

F 

effect  of  the  chance  failure  component. 


If  Eq.  (51)  is  considered  only  for  small  X  (this  implies  small 


\ ) ,  one  can  use  the  asymptotic  expansion  Eq.  (37)  for  P 
first  two  terms  only)  to  produce 


(take 


=  XT (2/u)/r (1/a) 


(l/\  -  ®) 


Using  either  Eq.  (51)  together  with  Table  I  or  Eq.  (52)  depend¬ 
ing  on  how  small  X  is,  the  value  of  k  ^  is  computed  and  plotted  in 

F 

Fig.  3  as  a  function  of  n  for  v  =  0.01  and  0.001. 


The  diagram  indicates  that,  with  all  other  parameters  fixed, k 

F 

increases  as  any  one  of  n,  y  and  a  increases.  Therefore,  the  decrease 

in  the  fatigue  sensitivity  k  ^  with  increasing  n  as  observed  above, 

C  /  \ 

gives  a  somewhat  fdlse  sense  of  security,  since  by  increasing  n,  k 

F 

increases.  For  example,  consider  the  case  of  >  =  0.01  with  a  =  3.  For 
(n )  *”2 

n  =  1.  k  1  '  <  10  and  therefore,  the  effect  of  the  existence  of  charce- 
F 

failure  component  may  be  neglected  at  least  for  the  estimation  of  ex- 

( n) 

pected  time  to  failure.  However,  for  n  =  100  and  1,000,  k  are  ap- 

F 

proximately  0.10  (n  =  100)  and  0.36  (n  =  1,000).  This  implies  that  the 
expected  time  to  first  failure  of  a  group  of  100  structures  is  10%  less 
than  that  of  a  group  of  100  (chance-failure-free)  reference  structures, 
and  36%  less  for  a  group  of  1,000  structures. 
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SECTION  IV 

DESIGN  AND  RELIABILITY 

The  relationship  between  the  conventional  design  procedure  and 
the  reliability  is  studied  in  this  section,  considering  an  (ideal¬ 
ized)  aircraft  wing  with  two  critical  locations  or  components.  It 
is  assumed  that  the  wing  root  on  the  upper  surface  (0^,)  is  subjected 
toa  possibility  of  failure  due  to  instability  (chance  failure)  and 
the  wing  root  on  the  lower  surface  (Op)  is  sensitive  to  fatigue  fail¬ 
ure. 


(a)  Gust  Data 


The  following  gust  data  in  the  flight  altitude  0  ~  10,000  ft. 
used  in  Ref.  4  are  employed  here  for  numerical  purposes.  The  dis¬ 
tribution  function  F(u)  of  the  gust  velocity  u  (ft. /sec.)  is  governed 
by 


F(u) 


.  -0.19u 

1  -  e 


(53) 


when  the  gust  is  observed  in  a  thunderstorm  turbulence  whereas 

,  ,  ,  ...  -0.344u  .  -0.208u 

F (u)  =  1  -  0 . 969e  -  0.031e 


(54) 


when  the  gust  is  observed  in  a  non-thunderstorm  turbulence  a.  -ociated 
with  general  operation. 


Furthermore,  it  is  observed  that  the  frequency  of  gust  occurrence 
is  14  gusts  or  7  gust  cycles  per  mile  in  a  thunderstorm  turbulence  and 
0.4  gusts  or  0.2  gust  cycles  per  mile  in  a  non-thunderstorm  turbulence 
(here  it  is  assumed  for  simplicity  that  an  up-gust  is  followed  by  a 
down-gust  and  vice-versa).  Therefore,  if  v  (u)  denotes  the  number  of 
gust  cycles  per  mile  with  velocity  larger  than  u,  then  for  a  thunder¬ 
storm  turbulence. 


v  (u) 


-0 . 19u 
7e 


o 

and  for  a  non-thunderstorm  turbulence. 


(55) 


V 

o 


(u) 


=  0.2  X  0.969e"°-344u  +  0.2 


x 


0, 031e 


-0.208u 


16 


(56) 


Under  further  observation  that  a  typical  aircraft  flies,  on 
the  average,  0.10%  of  the  total  flight  distance  within  a  thunder¬ 
storm  turbulence  while  18%  in  a  non- thunder storm  turbulence,  one 
can  combine  Eqs.  (55)  and  (56)  to  derive  the  number  of  gust  cycles 
per  hour,  v( u) ,  with  velocity  larger  than  u; 

v(u)  =  0.001V  x  7e  0.18V  x  0.2  x  0.969e 


+  0.18V  x  0.2  x  0.031e 


-0.208u 


where  V  is  the  flight  speed  of  the  aircraft. 

Assuming  V  =  400  miles/hour  in  the  present  study, 

.  .  _  -0.19u  _  „ ^  _  -0.208u  ,  .  -0.344u  , 

v  (u)  =  2.80e  +  0.432e  +  14e  (58) 

Since  v(u)  is  the  expected  number  (per  hour  of  flight)  of  gust 
cycles  with  velocity  larger  than  u  and  since  it  is  a  small  number 
compared  with  unity,  one  can  interpret  this  as  the  probability  that 
such  gust  cycles  will  occur  in  one  hour  of  flight. 

★ 

The  relative  frequency  v  (u)  of  gust  cycles  with  velocity  larger 
than  u  can  be  obtained  as 

*/  \  i  /  m\  n  -i  cn  “0.190u  -0.208u  -0.344u 

v  (u)  =  v(u)/v(0)  =  0.362e  +  0.025e  +  0.B13e 


The  values  of  v(u)  and  v*(u)  are  tabulated  in  Table  II  and  v*(u) 
is  plotted  in  Fig.  4  as  a  function  of  u. 

(b)  Ultimate  Strength  Design 


With  these  data  at  hand,  consider  the  design  of  a  wing  under  the 
specified  limit  gust  load  66  ft. /sec.  with  a  safety  factor  1.5. 


A  major  assumption  here  is  that  the  ultimate  failure  will  occur, 
due  to  upward  gust,  at  the  wing  root  on  the  upper  surface  in  the  form 
of  instability  rather  than  at  the  wing  root  on  the  lower  surface.  Such 
a  situation  can  arise  when  the  necessity  of  assuring  a  reasonable  fa¬ 
tigue  performance  results  in  a  design  of  the  wing  root  on  the  lower 
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surface  with  an  ultimate  strength  corresponding  to  a  safety  factor 
larger  than  1.5. 

Under  this  assumption,  the  chance  failure  of  the  wing  is  con¬ 
trolled  by  the  possible  instability  on  the  upper  surface. 

Hence,  the  limit  gust  velocity  66  ft. /sec.  and  the  nominal  ul¬ 
timate  strength  UFLN  (in  terms  of  the  acceleration  due  to  gravity  g) 
of  the  wing  is  related  by 

1. 5(1.0  +  66  Kn)  =  UFLN  (60) 

where  is  the  gust  sensitivity  factor,  a  function  of  flight  veloc¬ 
ity,  lift  coefficient  gradient,  effective  wing  area,  and,  the  number 
lt0  in  the  parentheses  represents  the  effect  of  the  level  flight. 

Suppose  that  the  wing  is  designed  sc  that  K  =  0.035  (a  repre- 

n 

sentative  value  indicated  in  Ref.  4)  and  UFLN  =  5.0  satisfying  Eq. 

(60)  . 

It  is  important  to  note  at  this  point  of  discussion  that  Jab- 

5  6 

lecki's  study  and  a  more  recent  study  by  Freudenthal  and  Wang  in¬ 
dicate  possible  unconservative  discrepancies  between  the  nominal  ul¬ 
timate  strength  and  the  actual  ultimate  strength.  From  these  studies, 
it  is  observed  that  50%  of  (structural)  specimens  sustained  approxi¬ 
mately  95%  of  the  nominal  strength  or  0.95  UFLN  is  the  median  value 
of  the  statistical  distribution  of  the  actual  ultimate  strength. There¬ 
fore,  taking  0.95UFLN  =  4.75  as  the  actual  ultimate  strength  UFL,  the 
gust  velocity  u^  at  which  the  ultimate  failure  will  occur  is  computed; 

1.0  +  i^K  =4.75  or  uD  =  107  ft. /sec.  (61) 

From  Fig.  4,  one  can  evaluate  the  relative  frequency  vD*  of  gust 
cycles  with  velocity  larger  than  uD  =  107  ft. /sec.;  v^*  =  2.52  x  10 
Multiplying  this  value  by  v (0)  =  17.2/hour,  the  frequency  of  gust  cycles 
with  velocity  larger  than  u^  or  the  probability  of  ultimate  failure  is 
obtained  as 

p  =4.33  x  10  ^/hour  (62) 
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and  therefore 


Vc  =  1/(4. 3 3  x  10-9)  =  2.31  x  IQ8  hours  (63) 

In  passing,  it  is  pointed  out  that  this  value  of  V  agrees 
7  c  o 

with  V  =  7.14  x  10  hours  for  L-188  Electra  and  V  =  2.38  x  10 

C  c  7 

hours  for  L-749  Constellation,  estimated  independently  by  Lockheed. 

(c)  Fatigue  Design 

In  order  to  consider  a  numerical  example  of  fatigue  design,  a 
typical  S-N  diagram  for  fatigue  failure  of  Mustang  wings  with  20% 
of  the  ultimate  load  as  the  mean  load  is  taken  from  Ref.  8  and  re¬ 
plotted  in  Fig.  5.  The  wing  under  consideration  is  assumed  to  have 
this  particular  S-N  diagram.  The  alternate  load  S  about  the  mean 
load  is  also  given  in  terms  of  percentage  of  the  ultimate  load  and 
N  is  the  number  of  cycles  of  S  the  wing  can  sustain  before  a  fatigue 
failure.  If  one  justifiably  interprets  the  ultimate  load  used  in 
Ref.  8  as  the  nominal  ultimate  load  of  the  wing,  the  mean  load  equal 
to  20%  of  the  ultimate  load  is  consistent  with  the  result  of  the  pre¬ 
ceding  ultimate  strength  design  in  which  the  nominal  ultimate  load  is 
5(g)  whereas  the  load  due  to  the  level  flight  is  1(g). 

The  statistical  distribution  of  N  at  a  constant  load  level  S  is 
assumed  to  be  of  the  Weibull  type. 

V*)  =  1  -  exp[-  {<x-NoW^-No)f]  (64) 

where  N  ,  N  and  a  are  constant.  In  this  connection,  it  is  pointed 
o  e  ‘ 

out  that  the  S-N  diagram  in  Fig.  5  is  interpreted  as  the  S-EN  diagram 
with  3N  representing  the  expected  value  of  N,  although  this  is  not 
clearly  indicated  in  Ref.  8.  Moreover  it  is  assumed,  as  in  Ref.  1, 
on  the  basis  of  the  experimental  evidence^8  that  (1)  the  distribution 
function  of  the  fatigue  life  N  under  random  load  is  also  of  the  Weibull 
type,  (2)  the  Palmgren-Miner  rule  or  its  modified  version  can  be  ap¬ 
plied  to  estimate  the  expected  value  of  the  life  under  random  fatigue 
using  S-EN  diagram  and  (3)  for  mathematical  expedience, the  scale  parame¬ 
ter  a  of  the  distribution  function  is  assumed  to  be  a  constant.  For 
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the  present  study,  a  =  3  is  assumed.  This  is  equivalent  to  the 
assumption  that  the  standard  deviation  5  of  the  (common)  logarithm 
of  the  life  N  is^ 

6  =  EUog1()N}  =  n/(2. 3030^/3)  =  0.186  (65) 

It  is  to  be  pointed  out  that,  practically  for  all  cases  (whether 
it  is  constant  load  fatigue  or  random  fatigue)  the  standard  deviation 
6  is  less  than  0.186.  It  can  be  shown  that  this  fact,  combined  with 
the  further  assumption  that  Nq  in  Eq.  (64)  is  zero,  will  produce  a 
conservative  reliability  estimate  at  least  within  the  reliability  range 
of  interest.  The  assumption  Nq=  0  reduces  the  form  of  the  Weibull 
distribution  Eq.  (64)  into  that  consistent  with  the  failure  rate  given 
in  Eq.  (2)  with  conversion  of  the  number  of  cycles  (N)  into  the  time 
(t),  which  in  the  present  case  takes  the  following  form. 

N  =  17. 2t,  N  =  17.2V  (t,  V  in  hours)  etc.  (66) 

6  F  F 

Observing  the  S~N  curve  in  Fig.  5,  one  can  assume  that  the  al¬ 
ternate  load  S  less  than  (>%  of  UFL  has  no  contribution  to  fatigue 
damage.  Hence,  constructing  six  (6)  intervals  of  u  beginning  at  u  = 

6  ft. /sec.,  with  interval  length  of  12  ft. /sec.  (except  for  the  last) 
as  shown  in  column  (a)  of  Table  III,  the  relative  frequency  p^  of  the 
gust  in  each  interval  is  computed  using  Fig.  4  and  listed  in  column 
(d) .  The  center  of  each  of  these  intervals  is  given  in  terms  of  ft./ 
sec.  in  column  (b)  and  in  terms  of  %  UFL  in  column  (c) .  In  this  con¬ 
version  from  (b)  to  (c) ,  the  fact  is  used  that  UFL  =  5(g)  and  K  = 

n 

0.035  (g/ft ./sec . ) . 

To  estimate  the  expected  fatigue  life  EN  under  random  load,  the 

R 

Palmgren-Miner  rule  is  employed.  For  this  purpose,  the  expected  fa¬ 
tigue  life  EN.  associated  with  the  alternative  load  S.  indicated  in 
l  i 

column  (c)  is  found  by  making  use  of  the  S-N  diagram  in  Fig.  5  and 
listed  in  column  (e) . 

According  to  the  Palmgren-Miner  rule, 

m 

Y  p.EN  /EN.  =  k  (67) 

/  1  K  1 

i=l 
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where  m  is  the  number  of  load  levels  or  the  intervals  considered 
(m  =  6  in  this  case)  and  k  a  parameter  depending  on  material,  struc¬ 
tural  configuration,  etc. 

With  the  aid  of  Table  III 
5 

V  p./£N.  =  2 . ]  L,:10~7+  2.58x10~7+  7.39x10_8, 

i-1  11  (68) 

+  1.60x10"8+  5.06x10“9  =  5.64x10  7 

in  which  the  effect  of  the  last  interval  is  neglected. 

Assuming  k  =  1.0  in  the  present  study  (the  standard  Palmgren- 
Miner  rule) , 

ENr  =  1.77xl06  (69) 

which  by  virtue  of  Eq.  (66) ,  can  be  converted  into  the  expected  life 

ET  in  terms  of  hours 
F 

ET  =  EN  /17 . 2  =  1.03xl05  (hrs)  (70) 

F  R 

Making  use  of  the  relation  between  ET  and  V  given  in  Eq.  (44) 

F  F 

with  a  =  3.0, 

VF  =  +  L/a)  =  i-iSxlO5  (hrs)  (71) 

Finally,  it  follows  from  Eqs .  (63)  and  (71)  that  the  parameter 

v  =  V  /v  defined  in  Section  II  is 
F  C 

y  =  V  /V  =  1. 15x10 5/^.  31X108)  =  4,98xl0_4  (72) 

F  C 

At  this  point,  it  is  emphasized  that,  although  for  design  pur¬ 
poses,  the  use  of  the  Palmgren-Miner  rule  and  the  like  seems  unavoid¬ 
able,  the  validity  of  such  a  rule  must  be  carefully  examined  by  speci¬ 
men  and  component  fatigue  test  under  random  as  well  as  constant  load¬ 
ing  conditions. 

(d)  Structural  Reliability 

Certain  reliability  statement  can  now  be  made  on  the  basis  of 
the  preceding  result. 
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First  of  all,  the  reliability  L^n^(t)  for  a  group  of  n  such 
wings  is  given  by  Eq.  (21) .  Explicitly,  it  is  written  as 

LU)(t)=  Lc(n)  (t)LF(n)  (t)  =  exp[-n(t/Vc)-  n(t/VF)a]  (73) 

c.  8 

with  V  =  1.15x10  h'-'urs,  V  =  2.31x10  nours  and  a  =  3.0,  This  is 

plotted  together  with  L '  ' (t)  and  L  'n  (t)  in  Fig.  6.  It  is  inter- 

C  F  ,  ^ 

esting  to  note  that  the  straight  lines  for  '  (t)  (chained)  and 

L ^  (t)  (dashed)  are  inscribed  (from  below)  by  the  curves  for  L^°^(t) 
(solid) . 

Reading  from  Fig.  6,  the  time  to  first  failure  under  the  speci¬ 
fied  reliabilities  0.9,  0.99  and  0.999  are  listed  in  Table  IV. 

With  y  =  4.98x10  4  and  a  =  3,  X  =  >n^  and  Pc ^ (using  Eq. 
(37)) can  be  computed  for  various  values  of  n  (Table  V). 

Making  use  of  Eq.  (52),  the  chance  failure  sensitivity,  as  it 
is  defined  in  Eq.  (49) ,  can  be  shown  to  be  equal  to  X  multiplied  by 
F (2/3)/r (1/3)  =  0.493,  which  is  small  compared  with  unity.  This  im¬ 
plies  that  the  expected  time  to  first  failure  with  n  actual  wings  is 
practically  the  same  as  that  with  n  chance-failure-free  structures. 

It  should  be  emphasized,  however,  that  this  is  true  only  when 
the  expected  values  are  compared.  In  fact.  Fig.  6  clearly  indicates 
that,  if  the  group  size  n  increases  and/or  higher  reliability  levels 
are  desired,  possibility  of  chance  failure  becomes  a  more  dominant 
factor. 


Since,  obviously,  the  chance  failure  can  no  way  be  prevented  by 

inspection,  the  wings  have  to  be  inspected  for  fatigue  damage  on  the 

basis  of  the  following  (conditional)  reliability  L  *^(t)  involving 

(n)  ,  F 

the  mortality  distribution  f  (t)  in  Eq.  (25); 

F 


*(") 


(t)  =  1  - 


1 

ynn 


I  (na/V_)  (t/Vi?)a"1e"nJt/VC+(t/VF)  }dt  (74) 


F  o 

This  provides  the  probability  that  the  first  failure  "due  to  fa¬ 
tigue"  will  not  occur  before  t.  It  can  be  shown  that,  for  reliability 
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values  close  to  unity  (such  as  0.99,  0,999  etc),  the  conditional 
reliability  L  *  ^  (t)  is  practically  identical  with  L  'n^  (t)  (relia- 

•?  P 

bility  of  the  first  failure  in  a  group  of  n  chance-failure-free 
structures)  indicated  in  Fig.  6  (dashed  lines).  Hence  the  times 
to  first  failure  "due  to  fatigue"  associated  with  reliabilities  0.9, 
0.99  and  0.999  are  read  from  Fig.  6  and  listed  in  Table  IV  (in  the 
parentheses) .  These  values  may  be  interpreted  as  first  inspection- 
free  period  under  the  specified  reliabilities.  The  difference  be¬ 
tween  the  figures  without  and  within  parentheses  indicates  the  ef¬ 
fect  of  chance  failure.  As  pointed  out  above,  it  is  larger,  as  n 
increases  and/or  higher  reliability  levels  are  demanded. 

Other  informations,  particularly  some  of  the  expected  values 
can  also  be  obtained  without  difficulty.  Making  use  of  Eq.  (34)  and 
asymptotic  expressions  in  Eqs.  (39)  and  (41),  ET^,  E  T^and  E 
can  be  evaluated  for  various  values  of  n  (Table  V) .  This  result  in¬ 
dicates,  for  example,  that  the  conditional  expected  time  to  first 

failure  under  chance  mode  is  5.77x10  hours  for  n  =  1,000,  a  value 

4 

much  smaller  than  1.045x10  hours,  the  conditional  expected  time  to 

first  failure  due  to  fatigue.  However,  the  probability  of  having 

-2 

chance  failure,  4.45x10  ,  is  small  compared  with  that  of  fatigue 

failure,  1-4.45x10  ^ . 
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SECTION  V 


DISCUSSION 


(a)  In  the  preceding  sections,  an  idealized  structure  is  con¬ 
sidered  in  which  only  two  locations  have  possibilities  of  failure. 

In  reality,  however,  a  major  structure  or  structural  component  con¬ 
tains  a  number  of  critical  locations.  If  such  locations  and  their 
failure  modes  are  identified  using  the  engineering  judgement  based 
on  the  experience  as  well  as  on  available  theoretical  and  experi¬ 
mental  knowledge  (the  aircraft  industry  is,  in  fact  mostly  responsi¬ 
ble  for  this  phase  of  study  because  of  their  experience  and  famili¬ 
arity  with  the  structural  details) ,  the  failure  rate  of  the  structure 
may  be  estimated  as  the  sum  of  the  failure  rates  associated  with  in¬ 
dividual  failure  locations.  It  can  be  shown  from  the  addition  rule 
of  probability,  that  the  failure  rate  thus  obtained  is  a  conservative 
estimate.  However,  it  may  be  overly  conservative.  For  example,  if 
locations  0, ,  0, . ,0  are  subjected  to  possible  chance  failures 

with  probabilities  of  failure  p  ,p  . ...,p  per  hour  and  if  the  chance 

1  2  n 

failures  considered  are  due  to  the  same  sequence  of  load,  as  in  the 
preceding  sections,  these  events  are  completely  dependent  and  the  max¬ 
imum  value  of  p.  serves  as  the  rate  of  chance  failure  of  the  entire 

l 

structure.  This  argument,  however,  does  not  seem  to  apply  to  the  fa¬ 
tigue  failure,  since  the  experimental  evidence  indicates  that  the  de¬ 
pendence  between  fatigue  pr^  ;sses  of  two  specimens  is  weak  even  under 
the  same  sequence  of  load.  Therefore,  for  fatigue  failure,  it  is  high¬ 
ly  recommended  to  add  the  rates  of  fatigue  failure  <;t  possible  locations 
within  the  structure  to  obtain  a  conservative  estimate  of  the  rate  of 
failure. 


A  direct  consequence  of  this  discussion  is  as  follows;  (1)  the 

rate  of  failure  of  a  wing,  in  which  the  fatigue  failure  is  possible 

on  the  upper  surface  (in  addition  to  the  lower  surface)  due  to  the 

load  during  taxiing,  <;aking-off  and  landing,  is  given  by  l/vc+  u  x 

(t/V  )a-1/V  +  a*  (t/v  *)3  Vv  *  where  a*  and  V  *  are  the  scale  pa- 
F  F  F  F  F 
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rameter  and  the  characteristic  value  of  the  fatigue  process  of  the 
upoer  surface  (assumed  here  is  that  the  chance  failure  in  compres¬ 
sion  on  the  lower  surface  is  negligible)  (2)  the  rate  of  two  such 

wings  of  an  aircraft  is  1/V  +  2a(t/V  )  ^  /  V  +  2a*(t/V  )a  ^/V  * . 

C  F  Jt  r  F 

(b)  Also,  in  the  preceding  sections,  the  wing  root  on  the  lower 

surface  is  assumed  to  have  the  fatigue  failure  rate  of  the  form 

a(t/V  )U  1/V  which  is  zero  at  t  =  0.  This  is  not  a  correct  model 
F  F 

since  it  could  fail  at  t  =  0+  under  an  extremely  severe  load.  The 
probability  l/vc*  Per  hour  of  such  a  failure,  however,  is  assumed 
(Section  IV,  (b)  )  to  be  much  smaller  than  1/V  of  the  upper  sur¬ 
face  and  hence  on  the  basis  of  the  discussion  (a)  above,  it  is 
neglected  in  the  analysis. 

(c)  In  order  to  take  into  consideration  the  fact  that  the  ultimate 
strength  has  a  statistical  distribution,  95%  of  nominal  ultimate 
strength,  a  median  value  of  its  empirical  distribution,  is  used  in 
the  analysis  as  deterministic  strength  (Section  IV  (bj  ).  Evidently, 
this  is  an  approximation.  The  method  is  well  known  to  deal  with 
such  statistical  variation  of  the  resisting  strength  for  a  more  rig¬ 
orous  reliability  analysis  (for  example.  Ref.  4) .  However,  such  a 
rigorous  treatment  is  not  performed  here  because  (1)  in  view  of  the 
relatively  larger  statistical  scatter  in  gust  load  distribution,  the 
approximation  seems  reasonable  and  (2)  it  may  make  the  analysis  ex¬ 
tremely  lengthy  and  cumbersome,  if  not  unmanageable,  and  hence  ob¬ 
scure  the  essential  purpose  of  the  present  study. 

(d)  r  he  same  reasoning  is  also  valid  for  other  simplifying  assump¬ 
tions  made  in  the  preceding  analysis.  For  example,  the  following 
possibilities  are  tacitly  disregarded;  (1)  the  effect  of  the  order 
of  application  of  different  load  levels  (this  is  always  neglected 
if  th  Palmgren-Miner  rule  is  used)  and  (2)  the  randomness  in  the 
number  of  gust  cycles  per  hour.  Also,  (3)  a  quasi-static  approach 
is  employed  here,  although  it  is  recognized  that  a  more  rigorous 
dynamic  analysis  treating  a  turbulent  gust  velocity  as  a  stochastic 
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process  (for  example.  Refs.  12  and  13)  is  possible  and  advisable 
to  check  the  validity  of  the  quasi-static  analysis.  However,  it 
is  pointed  out  that  such  a  dynamic  approach  has  its  own  difficul¬ 
ties.  For  example,  usually  the  turbulent  gust  velocity  and  the 
response  are  assumed  stationary  and  Gaussian.  This  is  a  question¬ 
able  assumption,  particularly  when  one  is  interested  in  a  chance 
failure  involving  extreme  values  of  the  gust  velocity.  Furthermore, 
a  chance  failure  implies  the  notorious  first  passage  time  problem 
when  the  load  is  treated  as  a  stochastic  process.  Recent  effort^^’^6 
to  deal  with  these  difficulties  are  noted  here.  It  is  hoped  that  the 
further  study  in  this  direction  will  make  it  possible  to  apply  the  dy¬ 
namic  analysis  to  realistic  reliability  problems 
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SECTION  VI 


CONCLUSIONS 

On  the  basis  of  assumed  additivity  of  the  risks  of  failure, 
the  distribution  function  of  the  time  to  first  failure  of  a  group 
of  nominally  identical  structures,  each  with  two  critical  loca¬ 
tions  respectively  subjected  to  chance  failure  and  fatigue  fail¬ 
ure,  is  derived  with  a  careful  discussion  of  the  assumptions  in¬ 
volved. 

Various  statistical  quantities  relevant  to  the  time  to  first 
failure,  such  as  (conditional)  expected  time  to  first  failure  due 
to  chance  failure  or  due  to  fatigue  failure  are  defined  and  their 
analytical  expressions  are  derived. 

A  definition  of  fatigue  sensitivity,  somewhat  different  from 
the  one  given  by  A.  M.  Freudenthal,  is  introduced  together  with  the 
definition  of  the  chance-failure  sensitivity. 

A  realistic  example  involving  aircraft  wings  designed  according 
to  the  current  method  of  design,  indicates  that  when  the  group  size 
n  is  increased  or  the  desired  level  of  reliability  is  raised,  the 
chance  failure  becomes  a  dominant  factor  to  be  considered.  Because 
of  this,  the  usefulness  of  fatigue  sensitivity  or  the  chance-failure 
sensitivity  is  limited,  if  defined  on  the  basis  of  expected  values  of 
the  time  to  first  failure.  It  is  shown,  however,  that  the  (condition¬ 
al)  distribution  of  the  time  to  first  failure  due  to  fatigue  is  almost 
identical  with  that  associated  with  chance-failure-free  structures  for 
reliability  values  close  to  unity.  Hence,  the  fatigue  damage  inspec¬ 
tion  may  be  planned  on  the  basis  of  the  distribution  function  of  the 
time  to  first  failure  associated  with  these  chance-failure-free  struc¬ 
tures  . 

Application  of  the  present  method  to  the  estimation  of  the  relia¬ 
bility  or  to  the  derivation  of  the  distribution  of  the  time  to  first 
failure  under  more  realistic  situations  involving  more  than  two  criti¬ 
cal  failure  locations,  is  also  discussed. 
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APPENDIX 


VARIANCE  OF  TIME  TO  FIRST  FAILURE 
The  second  moment  of  the  time  to  first  failure  T^n^  is  given 


E{T(n)}2  =  J  t2n{l/Vc+  a/VF(t/VF)a_1] 


.  -n [ t/V  +(t/V_)  1 
a-1,  c  F  " 


dt  (A-1) 


which,  after  integration  by  parts  and  setting  t  =  nt/V^,  becomes 


E(T^}2  =  2(Vc/n)2  J  Tie-11- 


(A-2) 


Since  it  can  be  shown  that 


P  <n)E  T(n) 

c  c 


(A-3) 


Eq.  (A— 2)  reduces  to 

E{T(n)}2  =  2(V  /n)  P  (n)E  T(n)  =  2ET(n)E  T(n) 

L  L  L  L 

Therefore,  the  variance  Var  of  T^  is 


(A-4) 


Var  T^n^  =  E{T^n^}2  -  {ET^}2 

=  ET^^  {2EcT'n^-  ET^} 


(A-5) 


and  the  coefficient  of  variation  v  is 


v  =  (2  EcT(n)/ET(n)-  1)1//2 


(A-6) 


Since  ET^  =  P  ^  V  /n  lias  been  computed  in  Section  II,  only 

/_  \  c  c 

E  Tv  needs  to  be  computed  (from  Eq.  (A-3)  on  an  IBM  7094).  With 

Vo  /  |.j\  /  j 

ETl  and  ECT '  computed,  v  can  be  obtained  from  Eq.  (A-6).  Table 
VI  lists  the  values  of  v. 

Recall  the  fact  that  the  coefficients  of  variation  of  T  ; 

and  T  ^  are  respectively  unity  and  [T  (l+2/a)/r2  (1+1/a) -1]  both 

r 
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independent  of  n.  The  1 

0'29  and  °'24  respectively  * ZT-T  ^  ^  ° ’53' 

4  and  5. 

variation  of  T(n>  areTt  lndl°ateS  that  “•»  coefficients 

between  those  of  T  <">  and  -  (n) 


r  and  T 
c  F 
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TABLE  I 


Probability  of  Chance  Failure  P^  ^ 


(a)  y  =  0.01 


2.0 

3.0 

4.0 

5.0 

1 

0.00881 

0.00888 

0.00902 

0.00914 

10 

0.0275 

0.0405 

0.0496 

0.0562 

100 

0.0838 

0.173 

0.247 

0.304 

1000 

0.237 

0.569 

0.779 

0.885 

(b)  y  = 

=  0.1 

o 

• 

CN 

3.0 

4.0 

5.0 

1 

0.0838 

0.0849 

0.0864 

0.0875 

10 

0.237 

0.332 

0.393 

0.435 

100 

0.545 

0.816 

0.921 

0.963 

1000 

0.865 

0.994 

0.9997 

0.999999 

(c)  v  = 

o 

_ 

2.0 

3.0 

4.0 

o 

in 

1 

.546 

.570 

.584 

.592 

10 

.865 

.958 

.984 

.993 

100 

.981 

.9994 

. 99998 

- 

1000 

.998 

- 

- 

TABLE  II 


Frequency  v(u)  per  hour  of  flight  and 
relative  frequency  v*(u)  of  gust  cycles 
with  velocity  larger  than  u. 


m 

v*  (u) 

0 

17.2 

1 

20 

8 . 42xl0-2 

4 . 89xl0-3 

40 

1 . 52xl0-3 

8.82X10"3 

60 

3.33xl0-5 

1 . 94xl0-6 

80 

7.31xl0-7 

4 . 2  5xl0-8 

100 

1. 62xl0-8 

9. 41xl0-10 

120 

3.59xl0-10 

2 . 08xl0-11 

Gust  Spectrum  and  Expected  Fatigue  Life 


TABLE  IV 

Flying  Hours  Under  Specified  Reliability  L 


SSiSSI 

1 

'  -  ‘  '  '  1 

10 

100 

1,000 

5 . 38x10 ^ 

2.52xl04 

1.17xl04 

5.02xl03 

.9 

(  "  ) 

(  "  ) 

(  ) 

(5.38xl03) 

2.40xl04 

l.lOxlO4 

4.80xl03 

1.59xl03 

.99 

(  "  ) 

(  •*  ) 

(5.02xl03) 

(2.32xl03) 

1.12xl04 

5.02xl03 

1. 59xl03 

2.12xl02 

.999 

<  ) 

(5.20x103) 

(2.40xl03) 

(l.llxlO3) 

TABLE  V 

Values  of  X,  P  ET(n),  EcT(n)  and  E-T 


10 


100 


1000 


4.98x10 


2.31x10 


1.07x10 


4.98x10 


pc‘"> 

ET*n*  (hr) 

4.45xl0_4 

3 .03xl05 

2 . 06xl0~3 

4. 76xl04 

9. 55xl0~3 

2.21xl04 

4 . 4  5xlO~2 

1.03xl03 

E  T(n)  (hr) 
F 


TABLE  VI 


Coefficient  of  Variation  v  of 

(n) 

Time  to  First  Failure  T 


(a)  v  =  0.01 


2.0 


0 . 530 


0.550 

0.604 


3.0 

4.0 

5.0 

0.360 

0.290  i 

0.240 

0.382 

0.310 

0.268 

0.444 

0.418 

0.417 

0.634 

0.737 

0.818 

10 

100 

1000 


2.0 

3.0 

i 

4.0  5.0 

0.551 

0.403 

0.331 

0.289 

0.604 

0.519 

0.496 

0.492 

0.725 

0.797 

0.871 

0.918 

0.893 

0.989 

0.999 

0.999 

10 

100 

1000 


2.0 

3.0 

4.0 

5.0 

0.725 

0.632 

0.596 

0.579 

0.893 

0.931 

0.959 

0.974 

0.982 

0.996 

0.997 

0.997 

0.998 

0.998 

0.998 

0.998 

Fig.  2  Fatigue  sensitivity  k  as  a  function  of  group  size 


ultimate  failure  load  (UFJJ)  . 


Fig,  6  Reliability  L  (t) ,  L_ '  (t)  and  L_ (t) . 
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